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*7h ! 1. Introduction 

a., 

The heavy ion community is on the search of quark-gluon plasma already for more 
^H. than 30 years. For this purpose large amount of money was spent. Thus it is timely 

to ask the question: what have we delivered in exchange for this support? Could 
we produce a piece of quark-gluon plasma? Have we determined the properties of 
this matter? Could we arrive to results which will be worthwhile to mention after 

■ a few years has passed? 

5— i ' In this lecture we try to shed some light to a few of these questions. 

a | 

2. Time evolution of heavy ion collision 

The different reaction models offer even qualitatively very dissimilar pictures for 
the flow chart of a heavy ion collision (see Fig.l). Some of the models assume that 
the first collision of the incoming nucleons produce directly most of the hadrons 
appearing in the final state. In contrast other models are based on the production 
of a collective intermediate state which expands and hadronize in a later stage only. 
Some of them use even the concept of thermal equilibrized intermediate state. The 
main goal of heavy ion research is to find out whether this collective state is formed 
in the heavy ion collision, and if it is formed then what are its properties. 
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TIME EVOLUTION OF HEAVY ION COLL. 



incident partons are coherent 
PRIMARY COLLISION BETWEEN PARTONS 

partons are incoherent \ 

final state 

$ 

| ? evolution of masless quark-gluon matter? 

^ 

PRODUCTION OF HOT MQ MATTER 

U 

| EVOLUTION OF MQ MATTER 

JJ- 

HADRONIZATION OF MQ MATTER 

evolution of hadron matter final state 

final state 

At very final state the resonances decay 

Fig. 1. The flowchart of the heavy ion collision. Notation: MQ = massive quark 
3. Qualitative considerations 

When we try to draw a picture of the heavy ion collisions, it is enlightening to recall 
two earlier ideas (perhaps in a somewhat subjective way). 

V.N. Gribov pointed out that in the nucleon-nucleus reaction at the first col- 
lision the incident nucleon looses its low momentum components, since their cross 
section are the largest (Ref. The remaining "wounded nucleon" (A. Bialas in 

Ref. [ |j) will interact with the next nucleon in a modified way. 

According to the above considerations we can distinguish the following type of 
collisions: 

a) soft collisions: the dress of both incoming nucleons is intact; 

b) semi- hard collisions: one of the incoming nucleons is "wounded", while the 
other incoming nucleon has an intact dressing; 

c) both incoming nucleons are wounded. 
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One more fact: in a recent calculation it was shown that the transverse mo- 
mentum distribution of partons produced in the collision of bombarding and target 
partons has an exponential shape [ |!|] . Fig. 2 displays the m t distributions of the 
produced massive quarks obtained in a perturbative QCD calculation, The spectral 
shape, which is exponential in a wide m t range is a result of the structure of the 
reaction matrix element. Thus the observation of an exponential spectrum is by no 
means an indication of any thermalization. 

With these considerations in mind we may phrase the following problems: 

• Is it true, that the strange quark-antiquark pairs (ss are produced thermally in 
the fireball? 

• Or are they produced in the primary collisions? 

• In the latter case the ss pairs are not a signature for quark-antiquark plasma. 

Nevertheless the free recombination of massive quarks and antiquarks is a sig- 
nature for the formation of deconfined quark-antiquark matter. 




Fig. 2. Invariant transverse momentum spectra of massive quarks produced in 
p + p reaction at y/s = 200 GeV. The incident parton masses are zero, the produced 
quarks are heavy with effective masses Mq = 300,400,500,1000 MeV. Figure is 
redrawn from Ref . [ || . 
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4. The early concept of the quark matter 

At the beginning for the description of hadron matter — quark-gluon plasma phase 
transition we used the following ideas: 

• at the beginning of heavy ion collision a quark-gluon plasma with massless quarks, 
antiquark and gluons is produced; 



• a first order phase transition occurs, quasistatically satisfying the Gibbs criteria; 

• homogeneous distributions and a few averaged particle properties (e.g. fi q , fi s , T) 
describe everything 




y 

Fig. 3. Top panel: Normalized rapidity spectra of hadrons produced in Pb+Pb 
collision at Elab = 158 AGeV. Data are from Refs.[||, [|. Bottom panel: particle 
ratios normalized to unity at the mid-rapidity. The lines indicate: a, S + /X~; b, 
S + /$; c, P - /tt-; d, pr/K~; e, ^/tt". 



The BRAHMS experiment at RHIC confirms the strong rapidity dependence, the 
proton to antiproton ratio was measured at different rapidities as: p/p = 0.64±0.03 
at y = and p/p = 0.41 ± 0.03 at y = 2 [ §. 

It was observed also, that the expansion of the hot region is very fast, excluding 
the possibility of a slow quasistatic phase transition. 

What is the way out from this contradictory situation? 

One has to introduce a new concept instead of assuming an ideal 
massless quark-gluon plasma. Such an attempt is done in the ALCOR 
model introducing the massive quark-antiquark matter. 
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5. The ALCOR (ALgebraic COalescence Rehadronization) 
model 

In the ALCOR model [ g, § we assume, that just before the hadronization the dense 
matter can be described as a mixture of dressed up, massive quarks and antiquarks. 
The effective mass of the gluons at this point is much larger than that of the quarks [ 
, consequently the gluon fission into quark-antiquark pairs is enhanced and massive 
gluons disappear from the mixture. During the hadronization the quark numbers, 
as well as the antiquark numbers, remain constant. 

At this point it is proper to insert a short explanation. At very high temper- 
ature, say at T > 3 * T c we may assume a non interacting quark-antiquark-gluon 
plasma with zero mass (current mass) constituents. However, as we approach the 
critical temperature, the interaction between the constituents will become stronger 
and stronger. This intensive interaction is expressed by the large effective mass of 
the constituents, which automatically reduce the pressure, faster then the energy 
density. 

5.1. Basic equations of the ALCOR hadronization model 

According to the coalescence assumption the hadron numbers are proportional to 
the number of quarks from which they consist. For baryons and antibaryons we 
assume 

P = C p b q b q b q qqq 

A = CAb q b q b s qqs 

5 = C5b q b s b s qs s 

fl = Cn b s b s b s sss 

Here the normalization coefficients, bi, are determined uniquely by the require- 
ment, that the number of the massive quarks do not change during the hadroniza- 
tion. Due to this principle we have 4 ( 6 , if u and d are treated separately) conserved 
quantities in contrast to the 2 parameters i^baryom Estrange of the thermal models 

This is the basic assumption for all quark counting methods The equations 
showing that the number of a given type of quarks in the final hadron population 
must be equal to the number of this type of quarks in the hadronizing quark matter 
are as follows: 

s = 3S1 + 2S + A + ~K + T] 
s = 3£1 + 2S + A + K + t] 
q = 3p + E + 2A + K + n 

q = 3p + E + 2A + K + n . (2) 



P = Cpbqbqbqqqq 
A = CAbqbqbsqqs 
5 = C-bqbsbs-qss 
fl = CTTb^b^b^sss 



7r = C n b q bq-qq 
K = C K b q bsqs 
K = C-^bqb s qs 
T] = C v b w b s ss (1) 
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5.2. The consequences of the simplified equations 

Applying the trivial assumption that the coalescence coefficient for a given type of 
particle and antiparticle is the same : C x = CW , we obtain the following antibaryon 
to baryon ratios: 
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(3) 



Here we introduce the notation 



D = b ^ qs = £ f 4 ) 
b^b^s K 

and thus with this we can can write eq. (^|) in the simple form 

|_„.Z | = D .I (5) 

A p a A ft a 

Such simple relations as shown in eqs. (||)-(||) were first obtained for the linear 
coalescence model ( 6, = 1 for all bi ) in Ref. [ [10| . 

5.3. Coalescence equations without isospin symmetry 

The effect of isospin asymmetry also can be taken into account in ALCOR by 
the proper quark combinatorics in the rehadronization. After relaxing the isospin 
symmetry assumption we treat the u and d quarks separately. 

The coalescence equations for baryonic and mesonic flavor clusters now read: 



p + 


= C p b u b u b,iuud 


7T+ 


— C n + b u b^ud 


n° 


= C n b u bdbdudd 


7T° 


— C n o (b u b^uu 


A° + E° 


= Ca b u bdbsuds 
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= C K - b^bsus 
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= C=o 6 U b s b s u s s 
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= C'k° bd bgds 




= C H - &d b s b s ds s 


K° 


= C^o b^b s d s 


ft" 


= Cnb s b s b s sss 


n 


= C' v bsb s ss 
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Here again the normalization coefficients, b u are determined uniquely by the 
quark number conservation demand: 

u = 2p+ + n +A + E° + 2 £+ + S° + tt+ + (tt° + ?/)/2 + if + 

u = 2p~ +n° + A° +^ + 2Z~ +Z° + + {n° +r ] ')/2 + K- 
d = 2n° +p+ +A° + S° + 2E" + ~- +n- + (tt° + r]')/2 + K° 

d = 2n°+p-+A + S + 2S + + 5 + + 7r + + (7r +r ? ')/2+^° 
s = 3 A - +2 ET +2 S° + A° + E° + £+ + £" + K~ +77 

s = 3n + + 2S + + 2S +A° + S° + S _ +E + + J ft: + ii' + + r7 (7) 

In eq. (fO) 7r is the number of directly produced pions. (Note that most of the 
observed pions are created in the decay of resonances.) 

5.4- Consequences of the isospin asymmetric equations 
The isospin asymmetric ALCOR gives the following baryon ratios: 
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(8) 



Let us introduce the notations: 

bdhds _ K° b u hus _ K+ B _ K°/K° 

b^hds If ' buhus K- ' D K+/K- W 

Using these factors B, C and D, antibaryon to baryon ratios can be written as 
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The correction factor C shows the effect of isospin. For a system with equal 
number of neutrons and protons C = 1. For other cases it can be obtained by 
inserting the experimental values into the defining equation, eq. (|To|), or calculating 
it from a model. In the ALCOR C = 1.0 due to isospin symmetry. 

Addendum: in the complete ALCOR calculation the resonances belonging to 
the two lowest baryon multiplcts and to the lowest two meson multiplets are also 
included. 

6. Hadron ratios at midrapidity at RHIC 

If we are interested of particle multiplicities in a restricted rapidity interval, we 
have to use in the equations described above the quark numbers belonging to this 
interval and then the hadron numbers obtained from the coalescence equations are 
also those belonging to this interval. 



7. Comparison with experimental data 

7.1. Parameter free predictions of ALCOR 

The set of relations eq. (|^) and eq. ( |To| ) are parameter free predictions. They agree 
surprisingly well with the experimental results at SPS and RHIC energies on the 
ratios p/p, A/A, S/H, SI /SI and K~/K + , as it was shown in Ref. [pH|. 

7.2. The full ALCOR calculation 

For the full calculation we have to define the three input parameters of ALCOR [ 
0,1: 

• number of new uu and dd pairs is obtained from the number of tt (or h ) ; 

• number of new strange ss pairs is obtained from the number of K~ (or K~ 

• the stopping factor (i.e. the percentage of incident quarks stopped into the inves- 
tigated rapidity interval) is obtained from the ratio K + /K~ . 

Now one can use the ALCOR model for Au+Au collision at y/s — 130 AGeV. 
We choose 211 newly produced uu and dd pairs, and 105.5 ss pairs (which means 
f s = 2N S s/ (N U u + N d 2) — 0.5) and assume that 5.8 % of the total nucleon numbers 
stopped into the central unit of rapidity (which means 9 protons and 14 neutrons) 
in the Au+Au collision. We can calculate the rapidity densities and hadronic ratios 
in the central rapidity region, where experimental data exist. Table 1 displays 
published experimental data and the appropriate calculated ALCOR results. The 
ALCOR parameters were fixed from the first 3 measured data. 
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ALCOR 


Experimental 


STAR 




model 


data 


reference 


h~ 


280 


280 ± 20 


[12] 


K-/-K- 


0.14 


0.15 ±0.01 


[13] 


K+/K- 


1.14 


1.14 ±0.06 


[13] 


W lv + 


0.63 


0.60 ±0.03 


[13] 


A/A 


0.73 


0.73 ±0.03 


[13] 




0.83 


0.82 ±0.08 


[13] 


E-/v- 


0.015 


0.013 ±0.01 


[14] 


<P/K*° 


0.38 


0.49 ±0.13 


[15] 


p~/h- 


0.087 


0.070 ±0.002 


[16] 



Table 1. Hadron production in Au±Au collision at y 7 ? = 130 AGcV from the 
ALCOR model and the experimental data from STAR [ [H| [0| [lj| |l§. The 
first 3 rows are the inputs for fit the 3 necessary ALCOR parameters. 
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Fig. 4. Calculated (open symbols) and measured (full symbols) hyperon ratios 
for SPS (WA97) and RHIC (STAR) experiments, see Ref. [0. 

8. Energy dependence of physical quantities 

Since the rapidity density of charged particles is increasing with increasing energy, 
we can investigate the energy dependence of the ALCOR parameters. From avail- 
able AGS, SPS and RHIC data (see Fig. 5), one can determine the 3 necessary 
parameters and calculate the results from the ALCOR model. Fig. 6. displays 
the tendency of the stopping into the mid-rapidity region and Fig. 7 shows the 
energy dependence of the number of newly produced light quark- antiquark pairs 
(dN u jil dy) and the strange ones indicated by f s in a logarithmic energy scale. 
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■ 71" - SPS, Pb+Pb • n" - RHIC, Au+Au, STAR-QM01 




2 5 10 20 50 100 200 

EXPERIMENTAL DATA in MID-RAPIDITY s' /2 (GeV) 

Fig. 5. The three input (measured) data for ALCOR as a function of energy. 
Data are from the AGS, SPS and RHIC experiments. 




Fig. 6. Energy dependence of the stopping parameter in the mid-rapidity 
region. 
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New light quark production at mid rapidity 



Strangeness ratio 




Fig. 7. The energy dependence of the ALCOR parameters: the rapidity density, 
dN uu /dy = dN dd /dy, of newly produced light quark-antiquark pairs (left panel); 
ratio of the newly produced strange quark-antiquark quarks to non-strange ones, 
f s (right panel). 



9. Charge fluctuation 



In a recent paper [ 18 interesting observation was made on the problem of charge 
fluctuation. In earlier papers it was suggested that the ratio 



D 



{ <5Q 2 > 
< N ch > 



(11) 



is D « 3 for hadron gas in equilibrium, while it is D ss 1 for the quark gluon plasma. 
Here < SQ 2 > is the average of the charge fluctuation and < N c h > is the average 
value of the charge multiplicity. The measured D value is near to 3 . That would 
indicate that in heavy ion reactions not a quark-gluon plasma, but a hadron gas 
was produced. However, in Ref. [ |l8), it was shown, that in the coalescence scenario 
(ALCOR) one can also expect a D « 3 value. Thus one can conclude that the 
measured charged fluctuation is also in agreement with the basic assumption of the 
ALCOR, i.e. during the hadronization the number of constituent quarks and the 
number of constituent antiquarks are conserved. 



10. Conclusion 

In this paper we have shown that the assumptions of the ALCOR model are in 
agreement with the experimental facts. These assumptions are: a) the matter just 
before the hadronization consists of a mixture of massive quarks and antiquarks; b) 
the number both of the quarks and the antiquarks are conserved during the sudden 
hadronization; c) the hadrons are formed by coalescence process. 
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We have investigated the energy dependence of particle ratios and thus found 
that the Wroblinski factor (/ s ) has a definite energy dependence. This indicates that 
the strange quark-antiquark pairs and the light quark-antiquark pairs are produced 
with different mechanisms having different energy dependence. 

We quoted Ref. [ [l8| , where it was shown that the charge fluctuations calculated 
from the massive quark matter is also in agreement with the experimental results. 

Finally we may conclude, that we have not seen the signature of an ideal quark- 
gluon plasma, but we found that in the heavy ion collisions a new type of 
matter, the massive quark-antiquark matter is produced. 

Note added in proof: The quark coalescence model was misunderstood many 
times and it induced surprisingly strong misinterpretations. This situation is well 
demonstrated in the following criticism [ |2(| : 

"The emphases in refs. [ [l(| [n] are in the multiplicative aspect of the prob- 
abilities of having quarks and antiquarks in the same region of the phase space in 
their formation of hadrons. That results in an undesirable feature of s and s imbal- 
ance in the linear version [ |l^| , which is not satisfactorily resolved in the nonlinear 
version [ |ll[ by the introduction of unknown factors." 

On the contrary, in Refs. [[n], |ll| particle ratios are calculated in a way where all 
unknown factors drop from the results. In the nonlinear version the bi coefficients are 
normalization factors ensuring the different particle conservation rules. To interpret 
these normalization factors as "the introduction of unknown factors" is equivalent 
to stating that conservation of e.g. electric charge is an unknown rule. 
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